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Abstract 

A two-step hybrid perturbation-Galerkin method is presented for the solution of a 
variety of differential equations type problems which involve a scalar parameter. The 
resulting (approximate) solution has the form of a sum where each term consists of 
the product of two functions. The first function is a function of the independent field 
variable(s) x, and the second is a function of the parameter A. In step one the functions 
of x are determined by forming a perturbation expansion in A. In step two the functions 
of A are determined through the use of the classical Bubnov-Galerkin method. The 
resulting hybrid method has the potential of overcoming some of the drawbacks of 
the perturbation and Bubnov-Galerkin methods applied separately, while combining 
some of the good features of each. In particular, the results can be useful well beyond 
the radius of convergence associated with the perturbation expansion. The hybrid 
method is applied with the aid of computer algebra to a simple two-point boundary 
value problem where the radius of convergence is finite and to a quantum eigenvalue 
problem where the radius of convergence is zero. For both problems the hybrid method 
apparently converges for an infinite range of the parameter A. The results obtained 
from the hybrid method are compared with approximate solutions obtained by other 
methods, and the applicability of the hybrid method to broader problem areas is 
discussed. 
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1. Introduction 


Perturbation solutions to differential equations type problems have been useful in a wide 
variety of applications; and it has long been realized how useful computer algebra can be in 
forming perturbation expansions. However, typically the algebra becomes more and more 
tedious as higher and higher order terms are computed; and frequently the computational 
effort rises so fast from term to term that even with computational assistance very few 
terms can be computed. Such luxuries as determining the radius of convergence of the 
expansion are rarely allowed, let alone situations where the expansion parameter can be 
modified to in effect increase the radius of convergence [1]. Thus for cases where higher 
order terms may have a significant effect it is important to make as much use of the 
information contained in the lower order terms as possible. The hybrid perturbation- 
Galerkin method described herein seems to greatly extend the power and usefulness of the 
perturbation method without adding significantly to the computational effort. 

The hybrid technique was apparently first studied by Ahmed K. Noor and collabora- 
tors. Their series of papers [8,9,10,12,13] combine the perturbation method, the Galerkin 
method [3], and the finite element method (or other discretization techniques) to attack 
a variety of structural mechanics problems. Their “reduced-basis method” allows some 
nonlinear problems requiring thousands of degrees of freedom after discretization to be 
computed using nonlinear systems of equations with only four to ten unknowns. Noor 
and collaborators have also applied the same general principles sans discretization to some 
thermal and structures problems [7,11]. We refer to the same general principle as the hy- 
brid perturbation-Galerkin technique because in using the technique without discretization 
we do not think in terms of having large bases to reduce. 

We believe the hybrid technique can be adapted to a wide variety of problem areas. 
Previous work by the present authors demonstrates applications to slender-body problems 
which result in singular perturbation expansions [4], In this paper the main problem con- 
cerns a well-known quantum eigenvalue problem, namely the perturbed one-dimensional 
harmonic oscillator. This problem in first order perturbation theory is treated or referred 
to in many different elementary quantum mechanics texts. Other work yet to be published 
involves boundary layer problems and problems in which the perturbation expansion is 
made at two or more values of the expansion parameter. 

In Section 2 we give a general description of the method. In Section 3 we analyze a 
simple two-point boundary value problem to illustrate the technique. Section 4 gives a 
treatment of the quantum anharmonic oscillator problem followed by concluding remarks 
in Section 5. 
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2. Description of the Method 

The method used in this study is a two-step hybrid analysis technique which is based 
upon a perturbation expansion technique followed by a Bubnov-Galerkin or variational 
technique. We suppose that we have a differential equation with independent variable x 
which involves a scalar parameter A and that we desire the solution to this equation for a 
particular value or range of values of A. 

In the first step a perturbation solution is developed with A as the expansion parameter. 
The perturbation solution takes the form of a sum over terms, where each term consists a 
perturbation function (a function of x) times a gauge function (a preassigned function of 
A). The expansion may be singular or regular. For a regular expansion the set of gauge 
functions consists simply of {1, A, A 1 , A 8 , . . .}. The drawbacks of the perturbation method 
are i) that frequently the radius of convergence in A is finite, ii) that even within the radius 
of convergence the rate of convergence may be slow and a large number of terms may be 
required to gain the desired accuracy, and Hi) that considerable computational effort may 
be expended computing higher order terms. 

In the second step we will keep the perturbation functions but replace the gauge func- 
tions by new amplitudes which depend on A. In the Bubnov-Galerkin or variational tech- 
nique one seeks an approximate solution in the form of a linear combination of specified 
(known) coordinate functions (functions of x) with unknown coefficients (or amplitudes) 
which are functions of A. This technique can work very well if a good set of coordinate 
functions is chosen. If a large set of coordinate functions is used, the computational ef- 
fort can be large since the matrices involved tend to be full rather than sparse (as in the 
finite element or finite difference methods). The essence of the hybrid approach is that 
the coordinate functions for the Galerkin step are chosen to be the perturbation functions 
computed in the first step. This provides a basis which clearly is relevant for small values 
of A. We demonstrate that the hybrid results can be dramatically better than the pertur- 
bation results and that a set of coordinate functions can be useful far outside the radius 
of convergence for the perturbation calculation in which they originate. 

3. A Simple Example 

To introduce the method we consider the following simple two-point boundary value 
problem: The differential equation, 


F — AV'+A = 0, 

is to hold for y in the range [-1, l] with boundary conditions 


(3.1) 
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(3.2) 


K(-1)=V(1)=0. 


The exact solution to this problem, which is 


V 


= y + 1 - 


2(i_ e Mv+i)) 

1_ C 2A 


(3.3) 


exhibits boundary layer behavior near y = 1 for large positive values of A and near y = - 1 
for large negative values of A. The center portion of the function approximates a straight 
line with slope equal to one in the limit of large absolute value of A. Exact solutions for 
a number of values of A are shown in Figure 1. A generalization of this problem has been 
studied by Ferguson [2], 

We look for an approximate solution of the form 


^ = £<b-(yK>(A). (3.4) 

j=i 

where in the first step the vy’s are computed from a (regular) perturbation expansion 
about A = 0 and in the second step the 5y in ’s will be computed using the Bubnov-Galerkin 
method. 


Step 1: 

To take step one first, the solution to (3.1) together with (3.2) is written in the form 

^ = !>,•(!/) + 0(A" +1 ). (3.5) 

;=i 

This is substituted into (3.1), the differentiation is distributed onto the various terms, 
the left hand side is expressed in the form of a power series in A, and the coefficient of 
each power of A is set to zero. The immediate result is a set of second order differential 
equations independent of the parameter A. These equations plus the imposition of the 
boundary conditions on each Vj serve to determine the t/y’s. Each t/y is a polynomial in 
y of degree jf + 1, and thus they collectively span an n-dimensional space. The first few 
terms are 


”l 

= y 2 - 

i, 


t> 2 

= y{y 2 

- l)/3, 


U3 

= (y 2 - 

•i)(y 2 - 

1)/12, 

t>4 

= v(y 2 

-i)(3y 5 

1 - 7)/l80, 

vs 

= (y 2 - 

• i)(y 4 - 

4y 2 + 3)/360. 
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It is traditional to stop the calculation at this point. If the computation has been done 
for sufficiently large values of n, an analysis using the ratio or root tests will show that the 
radius of convergence of the perturbation series is 7T. This may have been anticipated from 
the exact solution since the denominator is zero for A = ±miri for positive integer values 
of m. The convergence is limited to A in the range ( — tt, 7r) even though the singularity 
occurs on the imaginary A axis. The error of the perturbation solution may be measured 
by the £ 2 norm of the difference of the perturbation solution for fixed n and the exact 
solution (3.3) divided by the L 2 norm of the exact solution. In Figure 2 the logarithm 
of this perturbation solution relative error is plotted vs. log(A) for several different even 
values of n. As to be expected, the higher the value of n the more abruptly the error norm 
rises as the radius of convergence is reached. 

Step 2: 

In step two, the Vj for j = 1,2, ...,n serve as coefficient (or interpolation) functions 
for the Bubnov-Galerkin method. The approximation (3.4) with the t/ ; ’s known and the 
Sj t n s unknown are substituted into the left-hand side of (3.1). We would like, of course, 
for this quantity, the residual, to be zero for all values of y in the range [—1,1]. Since 
this is impossible, we settle for n conditions which serve to force the residual to be small. 
There are many different criteria which can serve this purpose, but we choose the Bubnov- 
Galerkin criterion which says that the residual is orthogonal to each of the coefficient 
functions Vj. The result is a linear system of n equations 

yi a »,;(A) 6j in — 6<(A), (* = 1>2, . . . , n) (3*7) 

j=i 

where 


«.-,;(*) = f ^iiy) [vj{y) - Xvj{y)]dy, 


bi = A J^Vi(y)dy. 
For n = 2, the 6 fi „ are given by 


(3.8) 


. 60 A . 60 A 2 . . 

^ = 60?A*’ = 6cT+"a* (39) 

The use of computer algebra allowed us to determine that <5 l>n for * = 1 through n and for 
any value of n has the form 


, e1Lv )/,j 
= r. 


(3.10) 
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where 


. (2 n - 2* + 1)!! 

(2< +l)!!i! n ,n — J+1) ' 


We have no proof of this formula, but it has been verified through n = 15. 

We measure the error of the hybrid approximation by computing the relative L 2 norm 
of the approximate solution (3.6) compared with the exact solution (3.3). In Figure 2, the 
log of the error of the hybrid solution is also plotted vs. A for a number of even values of 
n. We see that n = 2 gives a reasonable approximation for |A| less than about 0.5. Use of 
higher and higher values of n allows reasonable results for higher and higher values of |A|. 
In comparing the errors of the hybrid solutions with those of the perturbation solutions 
we see that for given values of n the hybrid solution are always more accurate. Further, 
while the radius of convergence has a dramatic effect of the perturbation solutions, there 
is no trace of an effect on the hybrid solutions. 

It was mentioned above that the denominator of the exact solution (3.1) has zeros for 
A = ±mni for positive integers m. It is interesting to note that these singularities are 
reflected in the 6j <n ' s. Since the 6 jn ’s are rational functions of A 2 , their singularities must 
correspond to the zeros of their denominators, the A„’s. The zeros of the A„’s, which are 
all purely imaginary, are plotted in Figure 3 for several low values of n. We see that jus n 
gets larger and larger the roots approach integer multiples of ni. 

For purposes of investigating convergence properties we have made extensive use of 
computer algebra. For less trivial problems, we may be forced to use numerical quadrature 
in computing the coefficients in the Galerkin equations, the a,/s and bj's introduced above. 
In that case it may be appropriate also to solve the system of Galerkin equations, similar to 
(3.7), numerically. In fact, for some problems it may be necessary to discretize the system 
of differential equations and to perform even the perturbation expansion numerically. Thus 
the mix of symbolic and numerical computation is problem dependent. However, the basic 
results seem to hold that i) the perturbation results can be dramatically improved by use 
of the second step in the hybrid method, and u) the radius of convergence limitations of 
the perturbation method do not extend to the hybrid method. We should add that for 
this problem the perturbation step has provided very little information. The perturbation 
functions Vj all contain a factor of y 2 — 1 due to the boundary conditions; and the set of 
functions 


w i = (»’ ~ 1) V* 


(3.12) 


spans the same function space as the t> } ’s. Moreover, the tz/ 3 ’s can be written upon in- 


5 


spection of the problem and without any computation. However, for many other problems 
the perturbation step plays a significant role in determining a small but appropriate set of 
coordinate functions to use as input for the Galerkin technique. In such problems a wild 
guess as to a set of basis functions might require a very large number of coefficients in the 
Galerkin equations to be determined in order to obtain the accuracy desired. 


4. Quantum One-Dimensional Anharmonic Oscillator 

The second problem is of a very different nature. The classical harmonic oscillator 
(mass on a weightless spring without damping) satisfies the equation 


mx — —k x. 


(4.1) 


Throughout time the sum E of the kinetic energy T — m x J /2 and the potential energy 
V 0 = kx 2 / 2 remains constant. 


mi 2 k x 2 

t + v q = — + — = e. 
2 2 


(4.2) 


Any small modification which makes the forcing term (the right hand side of (4.1)) non- 
linear results in an anharmonic oscillator. We then write the energy equation (4.2) as 


mx 2 kx 2 T , _ 
- + — + V = E. 
2 2 


(4.3) 


The quantum one-dimensional anharmonic oscillator satisfies the non-dimensionalized 
Schrodinger equation 


(H 0 + V[x))r/j(x) = +y + K(x)]v>(x) =Et/;(x) 


(4.4) 


with 


/ | yl>(x) | 2 dx = 1 . (4.5) 

J -oo 

We assume that x J /2 + V(x) goes to positive infinity as | x | becomes infinite. Then, the 
only solutions are bound states with 


V»(±oo) = 0, 


and E has a discrete spectrum. 

For purposes of our sample problem we choose the case 

V (x) = A x 4 with A > 0. 


(4.6) 


(4.7) 
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The solutions consist of energy eigenvalues E m and their corresponding wavefunctions 

(eigenfunctions) for m = 0, 1,2, For the harmonic oscillator, A = 0, the solutions 

are 


E m — 


= m+ i 


rpm{x) = <f>m(x) = hm{x)e lS / 2 , 
where the Hermite polynomial of degree m, h m (x), is given by 


(4.8) 


M*) = -7^ 

\J 2 m m! y 




(4.9) 


The (f> m alternate between being even and odd functions of x. It can easily be shown that 


[ <^m( x ) finix) dx — (4.10) 

J -oo 

where in this equation only, represents the Kronecker 6 function. 

For A 0, we believe that closed form solutions do not exist. Many elementary 
quantum mechanics textbooks discuss procedures for quantum perturbations and apply 
them to this problem for sake of illustration. Step one follows a standard procedure. 

Step 1: 

Let 



n- 1 

tM x ) = Y, A'>m,i( x ), 0m,O = <f>m, 


t=0 

with 

f <t>m{x) dx = 0 (t = l,2,...,n- 1) 

J — OO 

and let 

n-1 


^ Em t ty -^m t 0 = 

i=0 


Then substitute (4.11) and (4.13) into the differential equation 


(4.11) 

(4.12) 


(4.13) 


( 2dx 2 ^~ 2 "*’^ ac4 •®m)V’m( x ) — 0* (4-14) 

expand in powers of A, set the coefficients of A J to zero, and solve for E m j and rp m j(x). 

More specifically, the computation of the perturbation terms is as follows: The coeffi- 
cient of A 1 in the expansion is 
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( 2dx 2 2 — ( x + ■^m,l)0m,o(l) (4-15) 

Multiply both sides of (4.15) by <£ m (z), integrate over all z, use the normalization condition 
on <t> m {x) to get 


E m ,i= [ x 4 <f> 2 m {x) dx. (4.16) 

J -OO 

Now multiply both sides by <£„(: r) where n ^ m, integrate over all x, use the completeness 
conditions on the <£ m (z), and use (4.12) to get 


Wi = E (4. IT) 

r^in e "» e " 

To compute higher order terms in the perturbation expansion, note that setting the 
coefficient of A' in. the expansion of (4.14) to zero results in 




Id 2 x 2 


i-l 


2 dx 2 2 
which can be written as 


= (-x 4 + £ m ,i)t/\n,._i(z)+53 E m>i ^ mi , w (x) + E m>i ^ mi0 {x) (4.18) 

J =2 


( 2 dx 2 ^ 2 — fm,i{x) + E m ,i "4’m,o{ x ) (4-19) 

where fm, i is at this stage a known quantity given by 

t-i 

fm,i(x) = (~X + E mi x)tl> mi i-i(x ) + y: E m j i-j(y). (4.20) 

;=2 

Multiply both sides of (4.20) by 4> m {x) and integrate over all z. The left hand side vanishes 
and, consequently, 


/ oo roo 

4>m{x) dx — I <f> m (x) X 1 (z) dx. 

-oo J -oo 

The quantity rp m> i(x) evaluates to 


(4.21) 


V'm.i = 53 

n^m 


KM Uj 

€ n 


(4.22) 


The summations indicated in this procedure, while nominally over an infinite number 
of terms, are in practice all finite. For this reason, and since all the integrals cited above 
can be evaluated in closed form, the calculations can, in principle, be carried out to any 
order without any loss of accuracy. With the use of a computer algebra system such as 
MACSYMA [6] or Mathematica [14] it is relatively easy to get eight or ten terms. However, 
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by itself a high order perturbation computation is not very useful for this problem since 
the expansion has a radius of convergence in A which is zero! 

The fact that the radius of convergence is zero can be demonstrated by performing 
ratio and root tests on the series E m<i (t = 1,2,...). But this situation may have been 
anticipated by observing that for any value of A less than zero the solutions to (4.4) and 
(4.5) change drastically. We suddenly have all continuum states and no bound states. 

Before proceeding with step two it should be noted that the total wavefunction t/> m (x) as 
computed by this procedure has not yet been normalized (see (4.5)). The renormalization 
can be carried out as the last step in the perturbation calculation. However, it is not 
needed at all if we go on to step two of the hybrid method. 

Step 2: 

Despite the lack of convergence, we proceed to step two. For n terms in the approxi- 
mation let the wavefunction which is a perturbation of 4> m {x) be 

n-l 

= ) ] 6m,n,j( A) (4.23) 

1=0 

and let its associated energy eigenvalue be E m<n . The Galerkin orthogonality condition 
becomes 


r oo fid 2 x 2 - 

y_ Qo ^ ("2" = 0, (t = 0, 1, 


which may be written in the form 


n-l 


+ AC* m> i,/ — 0 (* — 0, 1, . . . , Tl — l) 

i= 0 


where 


/ OO 

xl>i(x)x ip j (x)dx 

-OO 

/ OO 

M x ) V»i(x) dx. 

-oo 

Nontrivial solutions are allowed only for 


,n- 1) 
( 4 . 24 ) 


(4.25) 


( 4 . 26 ) 


det(B m| n + A C mn E mn Dm,n) — 0 


( 4 . 27 ) 


w v 



where B, C and D are symmetric n-by-n matrices. Consider first the case of n = 1. Then 
the matrices B, C and D are all 1-by-l matrices. The quantity B evaluates to e m ; D 
is unity because of the normalization of the 4 > m ; and C evaluates to |(l + 2m + 2m 2 ), 


resulting in 



1 + 2m ( 3A(l + 2m + 2m 3 ) 
_ + - 


(4.28) 


For n = 2 the left hand side of (4.27) evaluates to a polynomial in m, A, and E m , n 
which is of total degree 2 in A and E m<n . There are two solutions, but the one which is of 
interest is the one which passes through E m<n = m + \ at A = 0. The form of Jj m>2 is 


o(m) + 6(m)A 

'c(m) + <f(m)A' 

2 . (c( m ) A ] J 

e(m) 

l e(m) J 

128e(m) 


(4.29) 


where a(m), b(m), c(m), d(m) and e(m) are all polynomials in m with integer coefficients. 
See the Appendix for further details. 

For higher values of n the left hand side of (4.27) evaluates to a polynomial of total 
degree n in A and E m , n • Thus (4.27) defines n curves in the A -E m ,n plane. However, the 
only curve of interest is the one which passes through the point (0, e m ). Examination 
of a number of specific cases indicates that the Taylor series of E m n as a function of A 
agrees with the first 2n terms of the perturbation solution. This is gratifying though 
not completely surprising, because it is known that the first n wavefunctions from the 
perturbation expansion can be used to determine an energy value which is correct to the 
first 2n terms in the energy expansion [5]. Thus the hybrid and perturbation methods give 
the same number of correct terms in the energy power series expansion. 

Once the energy eigenvalue E m>n has been determined, we may find the eigenvector 
6 mn (A) from which the wavefunction V» m n (x, A) may be constructed. As the last step the 
normalization is determined by requiring that 


f | VVn(*, A ) | J dx = 0. (4.30) 

J -OO 

While the Taylor series of the energy eigenvalues as a function of A may agree, there is 
a vast difference in the usefulness of the perturbation and hybrid results. The perturbation 
results, indeed, are virtually useless for the computation of energy eigenvalues and wave 
functions for any value of A, while the hybrid results seem to converge to the exact solutions 
as the number of terms increases. As in the previous problem we do not yet have a 
mathematical proof, but the empirical evidence in favor of convergence is strong. 

Consider first the energy eigenvalues. Figure 4 shows the results of the perturbation 
and hybrid approaches in computing the ground state (m = 0) energy level as a function 
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of A. The perturbation result for n = 2 is the same as the hybrid result for n = 1. Using 
higher perturbation terms seems to be counterproductive for almost all values of A on 
the scale shown, since at best only asymptotic convergence is exhibited. One the other 
hand, the hybrid results appear to be converging well for all A to the numerically obtained 
solutions. 

The numerical results were obtained by a shooting method and are presumed to be 
accurate to plotting accuracy. The initial point for the shooting method was a relatively 
large value of x with the requirement that the slope at the initial value of x be adjusted 
such that either the slope of (in the case of even m) or itself (in the case of odd m ) 
be zero at the origin. 

Figure 5 compares the hybrid results for the first four energy levels with numerical 
results. The errors always seem to diminish as the order is increased; and our work gives 
no hint of any limit to the range of convergence as A is increased. It may be noted, however, 
that the number of hybrid terms needed for a given level of accuracy does increase with 
energy. 

Convergence of the wavefunctions is considerably slower. Some representative cases are 
shown in Figures 6 and 7. 

While we are not able to predict how many hybrid terms are needed to gain a given 
level of accuracy for a given value of A, the convergence is monotone and the difference in 
the results for n — 1 terms and n terms seems to be a good indication. 

5. Concluding Remarks 

In this study we have treated two sample linear problems to exhibit some of the behavior 
of the hybrid method. The first is a simple two-point boundary value problem with a known 
exact solution. The perturbation method exhibits a radius of convergence equal to x. The 
hybrid method converges faster and appears to converge for all values of A. In common 
with the perturbation method, accurate solutions can be achieved with a small number of 
terms for small A. However, as the magnitude of A increases, more and more terms must 
be used to achieve a given level of accuracy. Use of computer algebra has allowed us to 
give an algebraic formula for the solution of the resulting Galerkin equations for any order. 
Thus, a proof of convergence may not be difficult to obtain. 

The second problem is an eigenvalue problem where the eigenvalues and eigenfunctions 
are perturbed. This problem has simple solutions for A = 0, but the perturbation method 
exhibits a zero radius of convergence. Empirical evidence suggests that the hybrid solutions 
converge for all positive values of A. Once again, accurate answers can be achieved with a 
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small number of terms for small A, and more and more terms are needed as A increases. 
A future publication will demonstrate that a numerical solution can be found valid in the 
limit of large A, and that this solution can be coupled with the expansion at A = 0 to give 
much improved convergence properties. Again computer algebra has allowed us to work 
without round-off error and thus to get more “pure” comparisons. 

Computer algebra has been used here primarily to investigate convergence properties. 
However, computer algebra can be expected to play a major role in more practical calcu- 
lations, particularly in the perturbation step. 

The improvements brought about by step two of the hybrid method are considerable, 
and yet the computational effort of step two typically is not as great as for step one. 
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Appendix 


For the second order hybrid solution for the quantum anharmonic oscillator problem 
of Section 4, we record the following results: 


6m, 0,0 = m + j, 

6m, 0,1 = 0 = 0) 

6m, i.i = ^g^(828 + 966m + 1031m 2 + 130m 3 + 65m 4 ), 
c m , o,o = f(l + 2m + 2m 2 ), 

Cm, 0,1 = C m>1)0 = -Hi±2l(21 + 17m + 17m 2 ), 

Cm, 1,1 = 5fj(3708 + 11838m + 16747m 2 + 9948m 3 + 5299m 4 + 390m 5 + 130m 6 ), 

dm, 0,0 = 1> 

dm, 0,1 = d m> i )0 = 0, 

156 + 422m + 487m 2 + 130m 3 + 65m 4 


(A.l) 


(A.2) 


(A.3) 


dm, 1,1 = 

It follows that E m> 2 has the form (4.29) with 
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a = 2(1 + 2m)(492 + 694m + 759m 2 + 130m 3 + 65m 4 ), 

6 = 3(1932 + 6286m + 9195m 2 + 5948m 3 + 3299m 4 + 390m 5 + 65m 6 ), 
c = 32(1 + 2m)(21 + 17m + 17m 2 ), 
d = 48(111 + 347m + 472m 2 + 250m s + 125m 4 ), 

e = 4(156 + 422m + 487m 2 + 130m 3 + 65m 4 ) . 

The Taylor series expansion of (4.29) about A = 0 is 

E m , = ^ + ^-^y + ±^ + om 

which is in agreement with (4.28) since 

a-c . 1 

— - = 6m, 0,0 = m + 2 


(AA) 


(A.5) 


(-*•«) 
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(.4.7) 


b-d_ 3, 

g — ^m.o.o — ~(l + 2m + 2m ). 
For m = 0 , 1 , 2 and 3, (4.29) evaluates to 
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FIGURE CAPTIONS 

1 

1. Solutions of the simple two-point boundary value problem for several values of the 
parameter A. 

2. Comparison of errors for the perturbation and hybrid techniques showing the depen- 
dence on the number of terms n and the parameter A. The radius of convergence, 
R c , of the perturbation expansion equals ir. 

3. Roots of the denominator of the (approximate) hybrid solution as dependent on the 
number of terms n. This shows that for large n that the singularities are located at 
±mxi for positive integer m. 

4. Comparison of hybrid results for the ground state (m = 0) energy with perturbation 
and numerical results. The number of terms used in the expansion is n. Hybrid 
results - solid lines; perturbation results - dashed lines; numerical (“exact”) results 
- circles. Perturbation and hybrid results are the same for n— 1. The perturbation 
results have zero radius of convergence. 

5. Comparison of the hybrid results with numerical results for the four lowest states. 
Hybrid results - solid lines; numerical (“exact”) results - circles. The accuracy 
becomes less as the energy level increases. 

6. Comparison of the ground state (m = 0) wavefunctions for A = 1 as computed by 
the perturbation method (dashed lines), the hybrid method (solid lines) and by a 
numerical shooting method (circles). Perturbation and hybrid results are the same 
for n = 1. 

7. Comparison of the wavefunctions for the four lowest energy states for A = 1/2 as 
computed by the hybrid method (various dashed and solid lines) and by a numerical 
shooting method (circles). 
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